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Summary 

We  have  collected  seismic,  geologic,  gravity  and  potential  field  data  to  develop  a 
model  of  the  North  Aftican  lithosphere  via  an  integrated  approach.  We  have  naodeled  the 
gravity  data  along  three  transects  across  North  Afiica:  1)  from  north  central  Egypt  across 
Libya  to  the  Hoggar  Uplift,  2)  ftom  central  Tunisia  across  northern  Algeria  to  the 
Atlantic  coast  of  Morocco,  3)  from  north  central  Algeria  across  the  Hoggar  Uplift  to 
northern  Chad.  These  transects  represent  ^  major  propagation  paths  tetween 
earthquakes  and  nuclear  explosions  and  seismograph  stations  located  within  North 
Afiica.  The  gravity  models  are  constrained  wherevCT  possible  by  geologic  and  other 
geophysical  information.  Crustal  thicknesses  we  have  obtained  for  the  Hoggar  Uplift 
(38  km)  agree  well  witii  preliminary  results  of  seismological  studies  by  other  research 
groups.  Our  modeling  results  also  suggest  that  the  manfie  anomaly  under  the  Hoggar 
Uplift  is  relatively  small  (200  km  wide,  less  than  20  krn  in  vertical  extent)  compared  to 

mantle  anomalies  associated  with  other  hot  spots  and  rift  zones. 

We  have  used  our  modeling  results  in  exjunction  with  other  geophysical  studies 
to  construct  a  map  of  crustal  thickness  for  Africa.  In  most  places  away  from  the 
continental  margin,  the  crust  in  Ncalh  Africa  is  35-40  km  thick.  Crustal  thinning  is 
associated  with  the  Sirt  Basin  and  the  West  and  Central  Afiican  rift  system.  Moderate 
crustal  thickening  (~41  km)  is  associated  with  the  Atlas  Mountains. 

In  addition  to  our  gravity  studies,  we  have  compiled  a  collection  of  seismograms 
ftom  Worldwide  Standardized  Seismognqih  Network  (WWSSN)  stations  locat^  at 
regional  distances  for  magnitude  ^  5.0  earthquakes  occurring  within  North  Africa 
(between  1963  and  1988)  and  their  associated  aftershocks.  We  have  selected  these 
earthquake  sequences  for  study  since  many  of  the  mainshocks  have  known  focal 
mechanisms  determined  from  teleseismic  waveform  riKxieling  or  first  motion  studies  of 
regional,  teleseismic  and/or  local  network  data.  Many  of  the  events  produced  numerous 
aftershocks  of  varying  magnitudes.  We  have  digitized  over  200  regixal  waveforms 
from  these  aftershock  sequences.  Using  the  lithospheric  models  obtained  firorn  the 
gravity  modeling  as  first  approximations  of  the  velocity  structure,  we  have  modeled  the 
regional  waveforms  of  these  earthquake  sequences.  We  have  also  collected  delay  time, 
travel  time,  source  parameter  and  magnitude  information  for  these  j^uences  for  our 
seismic  data  base.  The  entire  geophysical  and  geological  data  base  is  available  at  the  web 

site:  http://www.geo.utep.edu/nafr/nafr.html 
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Introduction 

North  Africa  has  been  the  object  of  cwiaderable  oil  exploration  and  a  fair  number 
of  papers  ccmtaining  geological  inteipretations  of  upper  crustal  structure  have  been 
published.  However,  very  few  papers  contain  basic  data  which  are  relevant  on  the  scale 
of  a  lithospheric  investigation.  Our  task  has  been  to  develop  a  model  of  No^  African 
lithosphere  via  an  integrated  analysis  of  seismic,  potential  field,  and  geologic  data.  In 
particuto,  we  have  ctwicentrated  on  constructing  detailed  2-D  models  from  known 

earthquake  source  regions  to  key  seismic  monitoring  stations  in  the  region.  An 

outgrowth  of  this  effort  is  a  data  base  of  geological  and  geophysical  information  which  is 
available  to  the  scientific  community  through  electronic  access. 

Although  data  are  sparse  in  many  regions,  it  is  clear  that  North  Africa  has 
e}q)erienced  a  conqtlex  evolutitm  (Figure  1).  The  structure  of  the  lithosphere  presently 
must  bear  the  signature  of  these  tectonic  events,  and  seismic  waves  propagating  through 
it  must  be  affect  The  Pan-African  orogeny  (750-500  Ma)  may  represent  the  time  at 
which  the  shield  areas  framed  by  terrane  accretion  (Black  et  aL  1994).  TOs  could  be 
rimilar  to  the  fwmatirai  of  the  China  region  in  the  Paleozoic  and  Mesozoic.  During  the 
Paleozoic,  the  Caledonian  and  Hercynian  orogenies  affected  the  regirai  ^  framed  a 
series  of  l^ins  and  uplifts  in  a  manner  similar  to  the  evolution  of  the  midcontinent 
region  of  the  U.S.  during  this  same  praiod.  In  the  Cretaceous,  Africa  was  subjected  to 
widespread  rifting  and  nearly  split  apart  (Fairhead,  1992).  During  tiie  Tertiary,  the 
Alpine  orogeny  affected  North  Africa,  most  notably  forming  the  Adas  Mountains. 
Volcanic  centers  such  as  Hoggar  and  Tibesti  fonned  on  topographic  uplifts  which  are 
believed  to  related  to  mantle  plumes  (Figure  1). 

Our  approach  to  constructing  a  lithospheric  model  and  data  base  for  North  Afri^ 
has  been  to  integrate  information  such  as  well  logs,  gravity  measurranents,  aerranagnetic 
surveys,  heat  flow  measurements,  and  geologic  maps  to  first  develop  lithosphraic 
models  along  key  prerfiles  in  North  Africa.  We  thrai  used  these  models  as  a  starting  point 
for  modeling  seismograms  at  regional  distances,  refining  the  models  as  the  seismic  data 
required.  In  following  subsections  we  describe  our  data  collection  process.  We  then 
dijyngR  the  results  of  our  gravity  modeling  and  analysis.  We  a^  present  prelumn^ 
results  of  regional  seismic  wavefram  modeling  and  a  crustal  thickness  map  for  Africa. 
The  appendices  describe  the  data  bases  we  have  collected  and  how  they  can  be  accessed 
electronically. 


Methods,  Assumptions  and  Procedures 
Geological  Data  Base 

We  have  digitized  the  locations  of  numerous  getdo^cal  features  within  the  North 
African  region  to  aid  in  our  interpretation  of  gravity  and  seismic  data.  Most  of  the 
digitized  data  comes  from  about  six  different  geologic/tectonic  maps  we  have  collected 
for  Nrath  Africa.  Other  features  have  been  digitized  from  figures  found  in  professional 
journals.  These  features  include:  surface  exposures  of  Precambrian  rocks  in  the  major 
uplifts,  depth  to  Paleozoic  and  Mesozoic  basement  within  the  Sirt  Basin  of  Libya, 
surface  faults  within  the  Sirt  Basin,  the  locatirai  of  Pan-African  suture  zemes  and 
Cretaceous  and  Cenozoic  age  rifts,  and  the  boundaries  of  West  and  Central  African  rift 
basins.  Maps  of  these  features  are  shown  in  Appendix  A,  along  with  details  on  how 
these  outlines  can  be  accessed  electronically. 
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Figure  1.  Generalized  tectonic  map  of  northern  Africa 


Gravity  Data  Base 

We  have  collected  gravity  data  frcm  a  number  of  sources.  First  we  begM  with 
the  unclassified  data  firtxn  die  global  data  base  maintained  by  the  Defen^  Mapping 
Agency.  For  these  data,  we  have  complete  informatioi  including  elevation  and  there  are 
over  36,000  data  points.  To  this  we  have  added  data  fix>m  Libya  which  were  con^iled 
by  one  of  our  graduate  students  (Suleiman,  1993).  There  data  (over  5000  points)  come 
from  a  variety  of  industrial  sources,  and  thus,  complete  information  is  laclong.  We  do 
have  latitude,  longitude,  and  Bouguer  anomaly  value  (density  -  2.67  gm/cc).  Details 
concerning  there  data  are  available  in  Appendix  B.  ^  ^ 

A  third  and  very  important  set  of  data  was  obtained  in  the  form  of  a  5*x5' 

Bouguer  gravity  anomaly  grid  fiom  GETECH,  which  is  associated  with  Iteeds 
University.  There  data  are  for  the  regirai  betweai  0*  and  37*N  and  12*W  and  40*E.  We 
are  only  free  to  use  there  data  in  the  form  of  maps  which  are  contoured  at  intervals  of  10 
mGal  or  greater.  Gravity  station  locations  were  also  provided  so  that  we  can  clearly  see 
the  coistraints  that  there  data  provide.  In  additicm,  the  staticMi  elevations  were  used  by 
GETECH  to  correct  errors  in  the  ETOP05  topographic  data  bare,  and  we  are  using  this 
commaeial  product  in  our  study.  The  ETOP05  data  bare  is  available  fiom  the  National 
Geophysical  Data  Center,  and  the  corrected  data  are  directly  available  from  GETECH. 

We  have  also  obtained  a  global  topographic  data  bare  fitan  GETECH  called  DMT5, 
which  has  5'x5'  resolution. 

All  three  data  rets  have  been  merged  to  prepuce  the  Bouguer  gravity  map  of 
Northern  Afiica  (Figure  2).  Details  on  the  availability  of  there  data  can  be  found  in 
Appendix  B. 

Seismic  Data  Base 

We  selected  a  series  of  magnitude  ^  5.0  earthquakes  and  their  associated 
aftershocks  occurring  between  1963  and  1988  (the  major  years  of  operation  of  the 
WWSSN)  for  study  (Figure  3,  Table  1).  Many  of  the  mainshocks  have  fwal 
mechanisms  detamined  from  telereismic  waveform  modeling  or  first  motitm  studies 
(Table  2).  Some  sequences  produced  numerous  aftwshocks,  and  we  collected  data  for 
all  events  above  magnitude  3.5.  The  data  we  collected  consisted  of:  1)  locatiOT,  origin 
time,  arrival  time,  delay  time,  magnitude,  and  first  motion  information  from  listings  of 
the  IntematicHial  Seisnwlogical  Centre  or  die  U.S.  Geological  Survey's  Earthquake  Data 
Report,  2)  published  focal  mechanisms,  focal  dq)th  and  other  source  parameter 
information  determined  fiiom  teleseismic  or  local  network  informatitMi,  and  3)  waveform 
data  celled  frcMn  microfiche  for  WWSSN  stations  located  at  regional  distances  fiom  the 
events  of  interest 

The  event  locations,  depths,  station  azimuths  and  distances,  P  and  S  arrival 
times,  P  and  S  delay  times,  and  magnitudes  were  entered  both  into  a  PC  based  version  of 

the  EXCEL  spreadsheet/data  bare  software  and  to  an  ORACLE  data  bare  that  is 

accessible  firmi  our  North  Afiica  web  site.  Exaii^)les  of  ^readsheets  and  other  details 
related  to  this  information  are  given  in  Appendix  C.  Figure  4  shows  an  e^^ple  of  how 
we  have  compared  P-wave  delay  times  versus  azimuth  for  events  within  different  regions 
of  North  Africa.  Figure  4a  shows  a  fairly  standard  means  of  depicting  the  delays  for  a 
miignitiiHp  4.5  aftershock  of  the  1969  northern  Red  Sea  earthquake  sequence.  Figure  4b 
shows  a  plane  parallel  coordinates  plot  of  the  same  information,  which  we  believe 
provides  a  ntorc  "usct  firiendly"  approach  to  viewing  delay  time  and  magnitude 
information.  More  exan^les  of  plane-parallel  coordinate  plots  are  found  in  Appendix  C. 

We  have  conduct^  an  extensive  search  of  the  literature  for  focal  mechanism 
information  and  believe  we  have  collected  all  available  focal  mechamsms  (Table  2). 
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stations  operating  during  the  period  of  interest. 


Table  1.  Seismic  Events  (1963-1988)  in  North  Africa  of  Interest  to  This  Study 


Date 

Time 

Magnitude 

Feb.  21,  1963 

1724 

5.0 

Gulf  of  Sidra 

May  11,  1963 

0110 

5.8 

Gulf  of  Sidra 

Oct.  20,  1963* 

1259 

5.6 

French  nuclear  test 

Jan.  1,  1965* 

2138 

5.0 

M'Sila,  Algeria 

Jan.  26,  1967 

1611 

5.1 

Tunisia 

July  13,  1967 

0210 

5.0 

Algeria 

Apr.  17,  1968* 

0912 

5.1 

Boudinar,  Morocco 

Mar.  24,  1969* 

1154 

5.2 

N.  Red  Sea,  Egypt 

Mar.  31,  1969 

0715 

6.0 

N.  Red  Sea,  Egypt 

Mar.  31,  1969 

2142 

5.0 

N.  Red  Sea,  Egypt 

Apr.  8,  1969* 

1031 

5.2 

N.  Red  Sea,  Egypt 

Apr.  16,  1969* 

0812 

5.0 

N.  Red  Sea,  Egypt 

Apr.  23,  1969* 

1337 

5.0 

N.  Red  Sea,  Egypt 

Dec.  1,  1970 

0102 

5.2 

Tunisia 

Jan.  12,  1970 

0815 

5.1 

N.  Red  Sea,  Egypt 

Sept.  2,  1972 

1453 

5.3 

Gulf  of  Sidra 

June  28,  1972 

0949 

5.5 

N.  Red  Sea,  Egypt 

Nov.  24,  1973 

1405 

5.1 

Algeria 

Nov.  24,  1973 

1522 

5.3 

Algeria 

April  29, 1974 

2004 

5.0 

Nile  Delta,  Egypt 

Sept.  4,  1974 

0629 

5.5 

Tripoli,  Libya 

Dec.  9,  1976 

2121 

5.2 

Tripoli,  Libya 

Jan.  1,  1977 

2046 

5.7 

A’in  Soltane,  Tunisia 

Dec.  9,  1978 

0713 

5.3 

W.  Egypt 

Oct.  10,  1980 

1539 

6.1 

El  Asnam,  Algeria 

Oct.  10,  1980+ 

1733 

5.2 

El  Asnam,  Algeria 

Oct.  10,  1980+ 

1908 

5.0 

El  Asnam,  Algeria 

Oct.  13,  1980+ 

0637 

5.0 

El  Asnam,  Algeria 

Oct.  14,  1980+ 

1734 

5.0 

El  Asnam,  Algeria 

Oct.  30,  1980+ 

2338 

5.1 

El  Asnam,  Algeria 

Nov.  8,  1980+ 

0754 

5.1 

El  Asnam,  Algeria 

Dec.  7,  1980 

1737 

5.1 

El  Asnam,  Algeria 

Dec.  12,  1980 

1332 

4.8 

El  Asnam,  Algeria 

Jan.  15,  1981 

0425 

5.1 

El  Asnam,  Algeria 

Feb.  1,  1981 

1320 

5.2 

El  Asnam,  Algeria 

Oct.  27,  1985+ 

1934 

5.7 

Constantine,  Algoia 

Mar.  26,  1988+ 

1207 

5.1 

Gulf  of  Sirt,  Libya 

*Seismograms  for  this  event  have  been  digitized 

+Seismograms  for  this  event  have  been  digitized  and  modeled.  See  Appendix  D. 
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Table  2.  Focal  Mechanisms  for  M^.O  North  African  Earthquakes  (1963-1988) 


Dal£ 


Feb.  21, 1963 
Jan.  1,  1965 

July  13,  1967 
Mar.  31, 1969 

Jan.  12,  1972 
June  28, 1972 
Apr.  29, 1974 
Sept.  4,  1974 
Oct.  10,  1980 
Oct.  30,  1980 
Nov.  8,  1980 


Mechanism  Data  Reference 

strike.diD.rake  Usfid 


243.75.90 
70,65,68 

167.30.90 
330,50,5 
325,50,-40 
8,86,-10 
350,80,10 
310,47,-17 
297,37,-141 
225,54,83 
42,40,72 
41,64,75 
50,30,97 


first  motion 
first  motion 
first  motion 
waveforms 
first  modem 
waveforms 
waveforms 
waveforms 
waveforms 
waveforms 
first  motion 
first  motion 
first  modem 


McKenzie  (1972) 
McKoizie  (1972) 
McKenzie  (1972) 
Ben-Menahem  et  al.  (1976) 
McKenzie  etal.  (1970) 
Pearce  (1980) 

Pearce  (1977) 
Ben-Menahem  et  al.  (1976) 
Westaway  (1990) 
Deschamps  et  al.  (1982) 
Yieleling  et  al.  (1989) 
Yielding  et  al.  (1989) 
Ouyed  et  al.  (1983) 
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Figure  4.  a)  Plot  of  P-wave  delay  time  (observed  minus  predicted  arrival  time  finom  listings  of  the  International  Seisnmlogical  Centre) 
versus  azimuth  from  event  to  station  for  an  aftershock  of  the  1969  northern  Red  Sea  earthquake  sequence,  b)  Plane  parallel 
coordinate  plot  of  same  information. 


All  regional  wavefam  data  that  were  available  have  been  collected.  We  have 
digitized  204  acceptable  seismograms  (Appendix  D).  Figure  5  shows  an  example  of  an 
M=5.0  aftershock  of  the  October  1980  El  Asnam  sequence  recorded  at  Toledo,  Spain 
(TOL)  located  at  625  km  from  the  event.  During  the  seccmd  year  of  the  research  we  have 
been  modeling  these  waveforms  to  validate  and  update  the  preliminary  lithospheric 
models  obtained  firom  the  gravity  data.  To  date,  we  have  modeled  88  waveforms 
recorded  at  WWSSN  seismic  stations  PTO,  TOL  and  MAL  (all  located  on  the  Iberian 
Peninsula)  from  aftershocks  of  the  1980  El  Asnam,  Algeria  earthquake  sequence.  We 
have  employed  reflectivity  techniques  (Randall  and  Taylor,  written  comm.,  1993)  in  this 
modeling  process. 


Results  and  Discussion 
Gravity  Modeling 

C^ce  we  had  merged  the  GETECH  data  into  our  existing  database,  we  selected 
three  profiles  across  the  region  to  model  (Figure  2).  These  profiles  represent  important 
propagation  paths  between  earthquakes  and  seismograph  stations  located  within  North 
Afiica.  Profile  1  extends  firom  central  Egypt,  the  site  of  several  recent  earthquakes  and 
seismograph  stations  (such  as  HLW),  across  central  Libya,  where  we  have  good 
geological  and  geophysical  control  on  crustal  structure  (e.g.  Suleiman,  1993;  Nyblade  et 
al.,  1996),  to  the  Hoggar  Uplift  of  southern  Algeria  (site  of  station  TAM,  Tamanrasset, 
and  several  French  nuclear  tests  in  the  early  1960’s).  Profile  2  extends  across  the  Atlas 
Mountains  frcMn  the  Atlantic  coast  of  Morocco  to  cwilral  Tunisia.  Seismograph  stations 
are  located  along  this  profile  in  Morocco,  Algeria  and  Tunisia,  and  earthquakes  occur 
throughout  the  region.  Profile  3  extends  from  north  central  Algeria  across  the  Atlas 
Mountains  to  the  Hoggar  uplift.  This  profile  is  along  the  propagation  path  between 
northern  Algerian  earthquakes  and  TAM.  A  preliminary  velocity  model  has  been 
obtained  by  McNamara  et  al.  (1996)  firom  regional  seismograms  whose  waves  propagate 

along  paths  that  parallel  Profile  3.  .  .  .  ^ 

The  density  naodel  we  have  obtained  for  Profile  1  is  shown  in  Figure  6.  The 
upper  part  of  our  nttodel  reflects  the  basin  structvre  of  Egypt,  libya  and  Algeria  obtained 
firom  published  geological  and  geophysical  studies  of  these  basins  (e.g.  Goudara,  1980, 
van  der  Meer  and  Cloetin^,  1W3;  Guirard  and  Maurin,  1992;  Exxon,  1985),  primarily 
related  to  their  petroleum  potential.  Since  our  knowledge  of  the  upper  crast  below  the 
sedimentary  basins  and  the  lower  crast  is  poor,  we  have  treated  these  units  as  having 
uniform  densities.  We  have  chosen  to  model  the  lower  crast  with  density  of  3.0  g/cc, 
upper  crast  with  density  2.70  g/cc,  and  mantle  with  density  3.30  g/cc.  These  densities 
were  chosen  as  average  continental  density  values  where  we  have  little  velocity  ccmsiraint 
on  the  crust.  We  varied  thicknesses  of  these  units  only  as  required  to  match  the  long 
wavelength  features  of  the  observed  gravity  data.  Note  that  we  have  not  atten^ted  to 
match  every  short  wavelengdi  feature  of  die  i^ofile,  since  our  knowledge  of  the  shallow 
geology  is  incwnplete.  Geologic  knowledge  in  these  regions  below  about  3  km  is 
roarse,  and  below  5  km  is  vague  or  completely  lacking.  We  have  used  information  from 
fiftailtvl  regicMial  studies,  where  available  (Hoggar,  Sirt  Basin),  and  continental-scaled 
geologic/tectcMiic  maps  and  resources  elsewhere.  The  shallow  geology  also  should 
influence  regional  waveform  inropagation  less  than  the  deeper  crustal  structure. 

Profile  1  shows  a  relatively  uniform  boundary  between  the  upper  and  lower 
crust  There  is  a  g^ual  shallowing  of  the  crast/mantle  boundary  between  the  Hoggar 
Uplift  (38  km)  and  the  Egyptian  coast  of  the  Mediterranean  (26  km)  (Sweeney,  1995). 

A  similar  eastward  shallowing  of  the  crast/mantle  boundary  is  observed  in  northern 
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Figure  5.  Seismograms  of  magnitude  5.0  aftershock  of  the  October  1980  El  Asnam,  Algeria,  earthquake  sequence.  Seismograms 
was  recorded  on  the  long-peri<^  vertical  instrument  operating  at  Toledo,  Spain  (TOL).  Top  trace  is  eat-west,  middel  trace  is  north- 
south,  bottom  trace  is  vertical. 
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Figure  6.  Interpretation  of  gravity  data  along  Profile  1  trending  West-East  from  Hoggar 
domal  uplift  (Algeria)  to  Meditenanean  coastline  (Egypt)  (total  length=2620  km). 
Above:  Observ^  and  calculated  gravity  values  for  each  station  along  profile.  Below: 
Crustal  model  used  to  generate  calculate  gravity  values.  Moho  discontinuity  shown  in 
bold.  Numbers  are  densities  for  the  bodies  (g/cc).  Upper  two  units  which  are  not 
labeled  are  Mesozoic  and  Paleozoic  sedimentary  rocks  with  densities  of  2.5  g/cc  and  2.6 
g/cc,  respectively. 
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Libya  (Suleiman,  1993).  Tlie  crustal  thickness  of  the  Hoggar  Uplift  is  in  good 
agreement  with  the  velocity  model  of  McNamara  et  al.  (1^6)  and  recent  receiver 
function  studies  at  TAM  (E.  Sandoval,  pers,  comm.,  1996). 

The  zone  between  approximately  280  -  400  km  shows  a  very  pow  match 
between  the  observed  and  calculated  gravity  values.  We  have  not  tried  to  model  the 
observed  anomaly  because  no  gravity  stad(His  exist  in  this  area  of  the  gravity  data  set. 
The  same  holds  true  between  600  -  820  km  and  1220  - 1350  km  along  Profile  1.  These 
areas  correspond  to  the  least  well  modeled  anomalies.  The  shallow  geologic  structure 
(<5  km)  in  die  Hoggar,  Murzuk  and  Sirt  Basins  is  very  well  constrained.  In  these 
regions  we  feel  that  the  local  anomaly  is  best  accounted  fOT  and  that  we  can  reasonably 
test  the  response  of  our  nxxlel  to  deeper  structure. 

Profile  1  (Figure  6)  intersects  Profile  3  OPigure  7)  at  a  distance  of 250  kra  Both 
profiles  have  a  total  crustal  thickness  of  38  km  at  this  point  and  all  layers  are  of  the  samp. 
thickness.  The  profiles  cross  in  the  northern  sectitxi  of  the  Illumeden  Basin  where  the 
thickness  of  Mesozoic  and  Cenozoic  sediments  are  well  constrained  (Tectonic  Map  of  the 
World,  1985).  Within  the  Sirt  Basin  and  the  Mediterranean  platfcxm  (2150-2620  km. 
Profile  1)  we  have  modeled  a  series  of  horsts  and  grabens.  The  horst  and  graben 
structure  resulting  from  Cretaceous  rifting  is  well  known  in  the  Sirt  Basin  (Goudarzi, 
1980;  Suleiman,  1993).  This  type  of  structure  appe^  to  occur  on  the  Mediterranean 
platform  as  well,  as  suggested  by  the  series  of  gravity  highs  and  lows  seen  in  diis 
region. 

Figures  7  and  8  show  interpretations  of  Profile  3.  The  shallow  upper  crustal 
stmcture  is  constrained  fiom  knowledge  of  basin  geometries  and  thicknesses  (e.g. 
Schurmann,  1974;  UNESCO,  1971;  Itesquer  et  al.,  1988).  The  Moho  depth  at  the 
northern  end  of  the  profile  is  constrained  by  the  results  of  the  European  Geotraverse 
(Ansorge  et  al.,  1992),  The  crust  appears  to  rapidly  thicken  beneaA  Ae  Adas 
Mountains,  then  gradually  thicken  to  the  south.  Recent  seismic  refraction  studies  (e.g. 
Wigger  et  al.,  1992)  have  modeled  the  Adas  Mountains  with  an  asymmetric  root 
extending  only  to  38  km  depth  in  the  north.  The  discrepancy  between  our  results  and  the 
refi:action  stupes  may  results  from  low  mande  velocities  (7.7  -  7.8  km/sec)  detected  in 
the  refraction  studies.  These  low  mande  velocities  suggest  diat  the  mande  density  may 
be  less  than  3.30  g/cc.  Again,  the  Hoggar  Uplift  can  be  modeled  with  a  38  km  Aick 
crust.  This  is  much  smallo*  than  observed  under  the  Ethiopian  and  Kenya  dcxnes. 

Figive  7  is  a  density  model  without  an  upwelling  mande  plume  breath  the 
Hoggar,  while  Figure  8  shows  the  same  density  model  with  an  upwelling  mande  plume 
add^  Note  that  there  are  zones  along  the  profile  with  poor  matches  between  observed 
and  calculated  anomalies.  The  zones  between  950  - 1050  km ,  1055  - 1200  km,  and 
1525  -  1775  km  contain  no  gravity  stations.  Thus,  we  have  not  tried  to  model  the 
observed  anomalies  in  these  areas.  Profile  3  intersects  both  Profile  1  in  the  south  and 
Profile  2  in  the  north  (Figure  2).  Where  the  profiles  intersect,  the  layer  thicknesses  have 
been  kept  constant.  Figure  8  shows  that  a  mande  plume  with  a  density  of  3.20  g/cc 
extending  horizontally  for  only  200  km  and  vertic^y  fixrni  38  U>  58  km  depth  is  all  that 
is  necessary  to  model  the  observed  gravity  anomaly  of  the  Hoggar. 

Between  1000  and  1900  km  lie  the  Sahara  Basins  which  include  the  TlliTi  and 
Mouydir  Basins  along  with  the  interbasin  divides.  Dq>th  to  the  Precambrian  bedrock  is 
well  constrained  within  this  area  of  Algeria  (Schurmann,  1974).  At  1900  km  lies  the 
South  Adas  fault,  separating  the  Adas  Mountains  to  the  north  ^m  the  Sahara  platform  to 
the  south.  The  structure  of  the  fault  is  not  well  studied.  It  is  believed  to  be  either  a  near¬ 
vertical  revCTse  fault  reactivated  in  the  Alpine  Orogeny  ot  a  con:q)lex  boundary  zone 
where  lower  angle  thin-skinned  thrust  faulting  penetrates  into  the  Saharan  domain 
(Outtani  et  al.,  1995).  The  former  view  is  mwe  widely  accepted  and  we  have  chosen  to 
model  the  fault  in  this  manner. 
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Figure  7.  Interpretation  of  gravity  data  along  Profile  3  trending  South-North  from  Air 
massif  (Niger)  to  Mediterranean  coasdine  (Algeria)  (total  lengdi=2240  km).  Above:  as 
in  Figure  6.  Below:  as  in  Figure  6,  with  no  mantie  upwelling  beneath  Hoggar  uplift. 
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Figure  8.  Interpretation  of  gravity  data  along  Profile  3  trending  South-North  from  Air 
massif  (Niger)  to  Mediterranean  coastline  (Algeria)  (total  lengA=2240  km).  Above:  as 
in  Figure  6.  Below:  as  in  Figure  6,  with  mantle  plume  beneath  the  Hoggar  uplift. 
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The  Adas  Maintains  arc  a  Tertiary  feature  which  cmtains  Precambrian  basement 
iq)lift  as  modeled  We  have  modeled  the  mountain  root  to  be  41  km  in  depth  to  match  the 
observed  gravity  low.  Previous  seismic  studies  have  modeled  the  crust  of  the  High  Adas 
Mountains  in  Morocco  to  be  38  km  thick.  Profile  3  crosses  the  Atlas  Mountains  in  north 
central  Algeria  where  the  range  is  about  2  km  lower  in  elevation.  The  apparcrit 
discrepancy  in  crustal  thickness  values  between  our  ^vity  models  and  the  seisimc 
studies  may  be  the  result  of  low  mantle  velocities  wWch  are  not  accounted  for  with  our 
mantle  density  of  3.3  g/cc.  We  may  choose  to  model  this  area  in  the  fu^  with  a  lower 
mantle  density  value  to  examine  the  effects.  North  of  the  Adas  Mountains  arc  duust 
nappes  emplaced  in  the  Tellian  Adas  during  the  Alpine  orogeny.  One  thrust  nappe  is 
modeled  in  Profile  3  between  2150  and  2240  km.  The  gravity  high  produced  by  the 
model  matches  the  observed  anomaly  well. 

Profile  2  (Figure  9)  extends  for  1840  km  in  an  east-west  directirai  across  the 
^proximate  center  of  the  Adas  Mountains.  The  profile  stretches  fiom  the  Adandc  coast 
of  Motocco  to  the  Mediterranean  shoreline  of  Tunisia.  Crustal  thickness  is  oidy  31  km 
along  the  west  coast  of  Motocco  near  the  Adandc,  but  quickly  thickens  to  41  km  beneath 
the  High  Atlas  of  Morocco  beginning  at  about  275  km  and  extending  to  a  distance  of 
1 125  km.  Here,  the  crust  begins  to  thin  as  the  Mediterranean  platform  is  crossed,  where 
Paleozoic  and  younger  sediments  reach  thicknesses  of  9  km  at  1840  km.  In  the  High 
Adas  Paleozoic  strata  outcrop  from  75-200  km  along  the  profile.  The  crust  thins  from  36 
km  at  1 150  km  to  34  km  at  1 450  km  across  this  part  of  the  Saharan  Adas.  Crustal 
thickness  remains  at  34  km  to  1750  km  along  the  profile,  after  which  it  quickly  thins  to  a 
minimiitn  of  28  km  entering  into  the  Mediterranean  platform.  Profile  3  is  intersected  at 
1 150  km  where  the  crust  is  36  km  thick  beneath  the  eastan  mar^  of  the  Saharan  Adas. 
No  gravity  stations  are  present  between  775  and  830  km,  thus  this  portion  of  the  profile 
was  not  modeled.  The  isolated  gravity  high  between  1450  and  1600  km  is  likely  due  to  a 
shallow  intrasive  body  that  has  not  been  documented  in  the  geological/geophysical 
literature  available  for  this  region. 

Gravity  Filtering 

Along  with  the  consonction  of  three  lithospheric  gravity  profiles,  we  have 
undertaken  wavelength  filtering  of  the  gravity  data.  In  gravity  data  the  longer 
wavelength  anomalies  arc  caused  by  deep^,  regional  structure,  while  the  shortCT 
wavelength  anomalies  are  caused  by  shallower,  localized  geologic  structure.  Figure  9 
and  10  show  the  gravity  data  filtered  to  retain  only  wavelengths  longer  than  1(X)  km 
(Figure  9)  and  200  km  (Figure  10).  In  Figure  9  the  distinct  anomalies  associated  with 
the  rifted  basins  of  the  Central  and  West  Afiican  Rift  System  (CWARS)  arc  still  evident 
Another  distinct  feature  on  the  figure  is  an  antMnaly  which  approximately  follows  the 
trace  of  the  Pan-African  sutures  where  the  WAC  and  Sah^  craton  are  joined.  In  Figure 
10,  sharp  anomalies  can  no  longer  be  associated  with  the  individual  aspects  of  the 
CWARS  system,  although  higha  gravity  values  still  exist  beneath  the  central  arm  of  the 
CWARS,  suggesting  that  the  mantie  is  likely  less  dense.  The  anomaly  associated  with 
the  Pan- Africa  suture  is  no  longCT  evident  and  the  gravity  lows  associated  with  the  main 
dranal  uplifts  are  subdued. 

Crustal  Thickness  Map 

Figure  1 1  shows  a  map  of  crastal  thickness  based  on  an  integrated  analysis  of  the 
limited  seismic  constraints  and  Bouguer  gravity  data.  Note  that  the  crust  is  generally  35- 
40  km  thick  away  from  the  continental  margins.  ThCTe  is  a  moderate  amount  of  mstal 
thickening  associated  with  the  Adas  Mountains.  Crustal  dunning  is  associated  with  the 
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Figure  9.  Lithospheric  model  for  Profile  2  (western  Morocco  to  northern  Algeria). 
Numbers  are  densities  in  g/cc.  Upper  two  units  are  as  in  Figure  6. 
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Figure  10.  Filtered  Bouguer  gravity  map  of  northern  Africa.  Wavelengths  >  100  km  passed. 
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Figure  11.  Filtered  Bouguer  gravity  map  of  northern  Africa.  Wavelengths  >  200  km  passed. 
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Figure  12.  Crastal  thickness  map  of  Africa.  Contour  interval  is  5  km. 
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Silt  Basin  and  the  West  and  Central  African  rift  system  and  the  adjacent  Canmoon  line. 
The  lack  of  a  major  crustal  anomaly  associated  with  the  dc»nal  upUfts  such  as  the  Hoggar 
is  surprising,  and  may  suggest  that  the  domal  uplifts  of  Hoggar,  Tibesti,  and  Darfur 
represent  rising,  less  dense  mande  more  than  crustal  thickening.  Locatitms  of  the  three 
gravity  profiles  are  shown  in  Figure  12.  Note  that  they  span  the  approximate  northern 
third  of  the  African  continent.  Qnstal  thickness  determinaticHi  in  Ais  area  is  our  main 
contribution  to  this  map.  The  thinnest  crast  of  the  African  continent  lies  along  the  Red 
Sea  rift  (about  20  km).  The  thickest  areas  of  North  African  crust  appear  to  be  in  the 
WAC  area  surrounding  the  Hoggar.  Since  none  of  our  profiles  cross  the  WAC,  we  can 
only  say  that  the  crusM  thickness  of  this  region  should  be  greater  than  40  km,  based  on 
our  modeling  across  the  Hoggar  re^rai  where  some  crustal  thinning  has  probably  taken 
place.  Gravity  profiles  2  and  3  indicate  that  the  Adas  Mountains  may  have  a  root 
extending  to  a  depth  of  41  km. 

Modeling  of  Regional  Seismic  Wav^orms 

Our  seismic  waveform  modeling  has  included  modeling  of  the  entire  regional 
wavetrain  from  P  through  to  S  waves.  We  have  collected  over  300  reccxds  and  digitized 
over  200  records  of  seismic  events  occurring  from  1963  to  1988  which  were  recorded  at 
12  WWSSN  seismic  stations  (Figure  3).  We  are  using  older  events  in  our  study  in  order 
to  focus  on  areas  widi  litde  recent  seismicity  (post- 1985)  and/or  to  gain  propagation 
paths  not  covered  by  more  recent  digital  staticms.  The  events  select  were  all  above 
magnitude  5.0  and  most  contained  numerous  aftershock  sequences.  Many  of  the 
mainshocks  have  focal  mechanisms  determined  frcKn  teleseismic  waveform  modeling  or 
first  motion  studies  (Table  2). 

In  order  to  ntodel  regicxial  waveforms  we  first  developed  a  preliminaiy  1-D 
velocity  model  (Figure  13)  based  upon  our  ^vity  modeling  studies  and  previous 
geophysical  studhes  (Seber  et  al.,  1996;  Saida  et  d.,  1996).  The  preliminary  1-D  model 
was  input  into  a  reflectivity  code  (Randall,  1989)  to  produce  syndetic  seismograms 
which  could  be  compared  to  the  obs^ed  seismograms.  Through  a  trial  and  error 
approach,  the  best  1-D  model  was  developed  to  create  synthetic  seismograms  similar  to 
the  observed  seismograms.  At  this  point  in  our  study  we  are  trying  to  get  the  best  match 
to  the  full  regional  waveform  from  our  1-D  models. 

Our  initial  modeling  results  are  for  88  observed  seismograms  at  WWSSN 
stations  TOL,  MAL  and  PTO  located  on  the  Iberian  Peninsula.  These  events  are  aU  part 
of  the  1980  El  Asnam,  Algeria  aftershock  sequence  (see  Appendix  D  for  waveforms). 
MAL  is  located  ~500  km  from  the  seismic  events,  while  TOL  and  PTO  are  located  650  to 
1000  km  for  the  epicenters.  Over  90%  of  the  synthetic  seismograms  required  a  shift  of 
less  than  10  seccmds  to  match  the  initial  portion  of  the  waveftxm.  Our  initial  model 
contained  four  layers  and  a  total  crustal  thickness  of  35  km  (Figure  13,  top).  For  station 
MAL  we  have  revised  this  1-D  nxxlel  to  have  4  layos  with  a  total  crustal  thickness  of  30 
km  (Figure  13,  bottom).  Models  for  stations  TOL  and  PTO  still  seem  to  be  best  modeled 
with  a  total  crustal  thickness  of  35  km.  We  have  been  able  to  match  the  early  part  of  the 
waveform  (P,PL)  well,  but  the  models  severely  mismatch  the  surface  wave  an^litudes. 
Much  of  this  problem  is  related  to  the  oceanic  Meditorranean  crust  blocking  surface 
wave  propagation.  We  have  also  noticed  large  amounts  of  energy  on  the  tangential 
ccxnponent  This  probably  reflects  backscattering  fix>m  structural  changes  along  the 
propagaticm  path  and/or  dipping  layers. 

We  will  continue  to  revise  Aese  1-D  velocity  models  and  begin  investigating 
other  propagation  paths  across  the  Mediterranean  and  North  Afiica.  We  will  use  events 
from  Ae  Hoggar  regicm  (including  ftench  nuclear  tests)  recorded  at  stations  bordering 
Ae  Red  Sea  to  mvestigate  the  structure  of  gravity  Profile  1.  We  will  also  investigate 
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Figure  13.  One-dimensional  velocity  models  for  TOL  and  PTO  (top)  and  MAL 
(bottom).  Bold  dashed  line  is  S-wave  velocity,  bold  solid  line  is  P-wave  velocity. 
Fainter  dashed  line  is  density. 


propagation  paths  near  gravity  profile  3  to  test  tiie  structure  modeled  there.  We  have 
records  from  some  events  in  northern  Libya  and  in  the  Mediterranean  Sea  north  of  Libya 
which  were  recorded  at  stations  on  the  Ibman  Peninsula  to  the  west  and  near  the  Red 
Sea  to  the  east.  We  plan  to  investigate  these  paths  for  2-D  structure  using  a  finite- 
difference  code  (Sandmeier,  1990)  after  we  have  determined  the  best  averaged  1-D 
structure  for  this  region. 


Conclusions 

Modeling  of  regional  sdsmic  phases  will  continue  throughout  the  next  year,  as  P. 
Dial,  who  is  currently  supported  as  a  teaching  assistant,  completes  the  modeling  as  part 
of  his  dissertation  research.  Geophysical  and  geologic^  data  bases  that  we  have 
collected  for  North  Africa  are  available  via  our  internet  web  site:  http://www.geo.utep. 
edu/nafr/nafr.html. 
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Appendix  A  -  Documentation  of  Geological  Data  Base 

Mjq)s  of  North  African  geological  features  that  we  have  digitized  are  available  at 
die  UTEP  Department  of  Geological  Sciences’  web  site  under  the  North  Africa  Project 
subheading  (http://www.geo.utep.edu/nafir/  nafir.html).  The  sources  of  data  used  for  the 
maps  we  have  generated  are  given  below,  along  with  information  on  the  scales  of  the 
original  maps.  Because  the  mt^s  on  the  web  site  are  in  color,  we  have  chosen  not  to 
reproduce  them  for  this  appendix. 

Figure  2  shows  the  outlines  of  exposed  Precambrian  rocks  m  the  Hoggar,  Tibesti 
and  Darfur  uplifts  (heavy,  dark  lines).  These  exposures  were  taken  from  UNESCO 
(1986)  (1:5,000,000). 

TTie  website  http://www.geo.utep.edu/nafr/nafrgeo.html  contains  a  color 
topographic  map  with  structural  features  shown  in  white.  Early  Cretaceous  basins  and 
faults  taken  from  Guirard  and  Maurin  (1992)  (-1:55,000,000).  Basement  faults  in 
north-central  Libya  are  frtmi  van  der  Meer  and  Qoetingh  (1993).  Structural  features  of 
central  Ubya  are  from  Goudarzi  (1980)  (-1: 1,500,000).  Structure  in  Algeria  and  Chad 
is  from  Dautria  and  Lesquer  (1989)  (1 : 10,500,000).  Early  Cretaceous  basins  in  the 
ncitheastem  portion  of  North  Africa  are  firtxn.Van  Houten  (1983)  (-1:1,500,000).  Lines 
in  western  portion  of  Neath  Africa  are  structures  shown  in  Kirton  and  Bourennani 
(1984)  (-1:3,000,000).  Well  locations  used  in  our  gravity  analysis  and  in  constructic>n 
of  the  stmctural  contours  are  shown  on  this  map  by  red  circles.  Majca-  geotectonic  units 
of  West  Africa  (Lesquer  et  al.,  1990)  are  also  shown.  Green  contours  on  the  map 
indicate  the  depth  to  the  Paleozoic  based  car  Goudarzi,  (1980;  -1:1,500,(XX)). 

The  website  htpp://www.geo.utep.edu/nafr/nafrgTOl.h^  shows  a  color 
topographic  map  with  the  location  of  our  gravity  profiles  in  white.  The  Central  and  West 
African  rift  system  is  shown  in  red.  Rift  basins  in  Niger,  Chad  and  the  central  African 
Republic  have  been  digitized  from  Genik  (1992)  (-1:2,500,000).  Precambrian 
exposures  of  the  Hoggar,  Tibesti  and  Dafur  uplifts  (UNESCO,  1986)  are  also  shown. 

We  have  also  digid^  a  Mtrfio  contour  map  from  Research  Group  fw 
Lithospheric  Structure  in  Tunisia  (1992)  (-1:600,000)  that  was  used  to  help  constrain 
our  gravity  ritodels. 
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Appendix  B  -  Documentation  of  Gravity  Data  Base 

Two  ^vity  data  bases  have  been  con:q)iled  for  public  use  in  North  Africa.  The 
first  of  these  in  the  Defense  Moping  Agency  ^MA)  data  base  that  includes  36,410 
^vity  points  distributed  over  North  Africa  north  of  10  degrees  north,  lliese  data 
include  principal  facts  and  are  in  the  following  format:  field  #1  -  latitude  in  decimal 
degrees,  field  #2  -  longitude  in  decimal  degrees,  field  #3  •  Bouguer  anomaly  in  milligals 
using  a  2.67  ^cc  reduction  density,  field  M  -  inner  zone  terrain  correction.  Hammer 
zxaies  (A-F),  in  milligals,  field  #5  -  outer  zone  terrain  correction.  Hammer  zones  G  and 
beyond,  in  milligals,  field  #6  -  elevation  above  sea  level  in  feet,  field  #7  -  observed 
gravity  in  milligals,  and  field  #8  -  source  numbo*  indicating  the  origin  of  the  data  point 
A  value  of  -99.99  milligals  in  the  terrain  correcticm  fields  means  that  the  terrain  correction 
for  that  point  has  not  been  made.  This  data  base  is  available  on  the  UTEP  web  page 
under  "DMA  gravity". 

The  second  data  base  available  is  mrxe  detailed  data  for  an  area  surrounding  the 
Gulf  of  Silt  Except  for  92  points  in  this  data  set  principal  facts  are  not  available  for 
these  data.  This  data  set  includes  5531  gravity  points.  The  format  is  the  same  as  for  the 
DMA  gravity  data  described  above,  excq>t  that  cmly  the  first  three  fields  of  data  are 
given.  This  data  base  is  available  on  the  UTEP  web  page  under  "Sirt  gravity". 

A  third  proprietary  data  base  was  used  to  prepare  the  maps  in  this  report  and 
those  available  on  the  UTEP  web  page.  This  data  base  contains  477,298  gravity  points, 
including  those  virtually  all  of  the  data  points  in  the  two  data  sets  mentioned  above  and 
has  been  gridded  at  5  minute  intervals.  Neitho'  the  data  points  nor  the  gridded  data  are 
available  to  tiie  public.  Individuals  with  specific  interest  in  these  data  should  contact  G. 
Randy  Keller  at  UTEP. 
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Appendix  C  •  Documentation  of  Seismic  Data  Base 

We  have  collected  earthquake  listings  fiom  the  bulletins  of  the  Intemadcmal 
Seismological  Centre  and  the  U.S.  Geological  Survey's  Earthquake  Data  Reports  for 
earthquakes  of  M^5.0  and  associated  aftershocks  occurring  between  1963  and  1988 
within  North  Afiica.  These  data  have  been  organized  into  a  data  base  using  the  PC  based 
program  EXCEL,  as  well  as  available  through  our  Internet  site  through  the  ORACLE 
data  base  management  program.  Diskettes  of  the  EXCEL  formatted  data  base  are 
available  upon  request. 

Figure  Cl  shows  an  example  of  the  main  spread  sheet  we  have  created  for  all 
events  of  interest  (example  in  EXCEL).  The  spread  sheet  contains  informaticm  on  the 
date,  time,  latitude  (lat),  longitude  (Ion),  depth,  magnitude  (Ms,  mb  or  other),  and 
location. 

For  each  earthquake  for  which  we  have  obtained  a  listing,  we  have  created  the 
^readsheet  (exan^le  in  EXC^)  shown  in  Fi^e  C2.  The  first  line  of  the  spreacisheet 
ctmtains  header  informaticm  on  the  earthquake  itself.  Then  we  list  all  stations  within  20* 
of  Ae  earthquake  and  note  the  station  distance,  statitm  time,  P-wave  arrival  time,  S-wave 
arrival  time,  and  arrival  time  of  otiier  phases,  if  given.  We  also  record  the  delay  time 
(observed  arrival  time  minus  predicted  time).  Magnitude  informaticm  (usually  log  A/T)  is 
dso  noted,  when  available. 

We  have  used  the  data  craitained  in  the  spreadsheet  data  base  to  generate  plots 
showing  travel  time  delay  information  fw  earthquakes  occurring  in  areas  of  Algeria  and 
Egypt  Using  plane-parallel  coordinate  projecticms  we  displayed  d^  base  information 
including  earthquake  parameters  such  as  distance,  azimuth,  travel  time  delay  and 
mflgnitiiHft  for  WWSSN  Stations  of  interest.  For  events  occurring  in  Egypt  we  used 
travel  tinv-  delays  fw  the  following  WWSSN  stations:  JER,  AAE,  HLW  and  1ST.  The 
stations  of  interest  for  Algerian  events  were:  LOR,  MAL,  TOL,  AGU,  TRI  and  PTO. 
Examples  of  these  plots  are  shown  in  Figure  C3  and  C4. 

In  addition  to  creating  these  plane-parallel  corstlinate  plots  we  have  also 
displayed  delay  informaticNi  in  map  view.  Examples  of  these  maps  are  shown  in  Figures 
C5  and  C6  for  HLW  and  PTO. 

Digitized  seismograms  are  shown  in  Figures  C7  to  C40.  For  each  figure  a  label 
is  shown  at  the  upper  ri^t  of  each  seismogram.  The  first  numbers  represent  the  start 
time  (hour  and  minute  UT)  of  the  seismogram.  This  time  is  followed  by  the  3  character 
WWSSN  station  code  and  a  1  character  component  type  (r=radial,  t=tangential, 
z=vertical).  Seismograms  are  available  upon  request  from  P.  Dial  (dial@geo.utep.edu) 
and  are  written  in  SAC  (seismic  analysis  code)  format. 
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Figure  Cl.  Example  of  page  from  master  event  spreadsheet  that  contains  information  on  all  earthquakes  we  plan  to  study. 
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Figure  C2:  Example  of  spreadsheet  data  for  an  aftershock  of  the  1969  northern  Red  Sea  earthquake  sequence. 


Figure  C3.  Plane-parallel  coordinate  projection  showing  earthquake  parameters  for 
Egyptian  events  recorded  by  WWSSN  stations  in  1969. 
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Figure  C4.  Plane-parallel  coordinate  projectiwi  showing  earthquake  parameter  variations 
for  WWSSN  stations  recording  Algerian  earthquakes  between  October  1980  and 
February  1981. 


Figure  C5.  Azimuthal  equal  area  projection  showing  travel  time  delays  for  the  WWSSN 
station  HLW  (triangle).  Circles  denote  negative  delays  and  crosses  indicate  positive 
delays.  Symbol  size  is  scaled  between  +/-  6.0  seconds. 
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Figure  C6.  Azimuthal  equal  area  projection  showing  travel  time  delays  for  Algerian 
earthquakes  recorded  at  PTO  (triangle).  Circles  and  crosses  are  as  in  Figure  C5. 
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Seismograms  of  El  Asnam  aftershocks,  October  10, 1980. 
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C8.  Seismograms  of  El  Asnam  aftershocks,  October  10, 1980. 


0540mal.r 


Figure  C9,  Seismograms  of  El  Asnam  aftershocks,  October  10, 1980. 
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Figure  CIO.  Seismograms  of  El  Asnam  aftershocks,  October  10, 1980. 
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Figure  C12.  Seismograms  of  El  Asnam  aftershocks,  October  10, 1980. 
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Figure  C14.  Seismograms  of  El  Asnam  aftershocks,  October  10, 1980. 
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Figure  C15.  Seismograms  of  El  Asnam  aftershocks,  October  10  (2127)  and  October  12  (0036),  1980. 


Seismograms  of  El  Asnam  aftershocks,  October  12, 1980. 
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Figure  C17.  Seismograms  of  El  Asnam  aftershocks,  October  12, 1980. 


Figure  C18.  Seismograms  of  El  Asnam  aftershocks,  October  13, 1980. 


Seismograms  of  El  Asnam  aftershocks,  October  13, 1980. 
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Figure  C20.  Seismograms  of  El  Asnam  aftershocks,  October  13, 1980. 


1434malz 


Seismograms  of  El  Asnam  aftershocks,  October  13, 1980. 
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Seismograms  of  El  Asnam  aftershocks,  October  13, 1980. 


Figure  C23.  Seismograms  of  El  Asnam  aftershocks,  October  13, 1980. 
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Figure  C24,  Seismograms  of  El  Asnam  aftershocks,  October  14, 1980. 
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Figure  C25.  Seismograms  of  El  Asnam  aftershocks,  October  14, 1980. 
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Figure  C26.  Seismograms  of  El  Asnam  aftershocks,  October  14, 1980. 


Os318mole.sac.O 


■■■V. 


56 


Figure  C27.  Seismograms  of  El  Asnam  aftershocks,  October  15, 1980. 
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Figure  C28.  Seismograms  of  El  Asnam  aftershocks,  October  19, 1980. 


Seismograms  of  El  Asnam  aftershocks,  October  21  (1651)  and  October  22  (1624),  1980. 
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Figure  C30. 
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Seismograms  of  El  Asnam  aftershocks,  October  23, 1980. 


Seismograms  of  El  Asnam  aftershocks,  October  23, 1980, 
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Figure  C33.  Seismograms  of  El  Asnam  aftershocks,  October  26, 1980. 
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Figure  C34.  Seismograms  of  El  Asnam  aftershocks,  October  30, 1980. 
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Figure  C35.  Seismograms  of  El  Asnam  aftershocks,  October  30, 1980, 
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Figure  C36.  Seismograms  of  El  Asnam  aftershocks,  October  30, 1980. 
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Figure  C37.  Seismograms  of  El  Asnam  aftershocks,  November  8, 1980. 
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Figure  C38.  Seismograms  of  El  Asnam  aftershocks,  November  8  (0209),  1980  and  of  Constantine  aftershocks,  October  27, 1985 


1211plo.r 


Figure  C40.  Seismograms  of  Gulf  of  Sirt  earthquake,  March  26, 1988. 
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